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ABSTRACT 


The communication channel between a satellite and a ship earth station (SES) 1s 
described by a model which includes multipath fading, doppler shift and noise. Multi- 
path fading is caused by reflections from the sea surface. These reflections can affect the 
system performance, especially at low elevation angles or when SES is using low gain 
antennas. Doppler shift 1s a very important effect when using low altitude satellites, be- 
cause of the high velocities involved. 

This thesis describes and presents a software simulator for multipath fading in the 
maritime communications environment. Analysis of throughput of an unslotted Aloha 


maritime mobile satellite communication channel is also presented. 
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I. INTRODUCTION 


Future satellite systems for maritime, aeronautical and land mobile applications 
(with mobile terminals for low data rate and digital voice transmission) will have small 
antennas having gains between 2 and 12 dB and a broad beamwidth. These antennas 
pick up multipath reflections of satellite signals causing severe fading of the received 
signals. 

Several fading simulators, using hardware, software, or both, have been designed. 
Among them are: JPL, CRC and Virginia Tech. All of them deal with land mobile 
communications and with geostationary satellites. 

In the future, low altitude satellites will become important alternatives for mobile 
communications because of the crowded geostationary orbit. Studies have demonstrated 
that five satellites in circular orbit will provide continuous single coverage of the earth; 
seven satellites will provide double coverage at some minimum elevation angle. By the 
same criterion it was found that eleven satellites can provide continuous triple coverage 
and fourteen satellites, quadruple coverage. Continuous single coverage means that any 
two stations on earth can communicate using this constellation of five satellites at anv 
moment. Double coverage, on the other hand, will provide an additional communication 
path; triple coverage will provide three different paths and so on. 

Besides the econonucal achievement of maximum coverage, many other factors have 
to be taken into account in system design. These results may prove of some interest to 
planners of satellite systems for mobile communications and surveillance [Ref. 1] 

Low earth orbits (LEO) are generally circular or slightly eccentric. In order to 
simplify the analysis and equations described in this thesis, the focus will be on circular 
orbits. 

This thesis presents a fading simulation program for mobile maritime satellite com- 
munications involving low altitude satellites. 

Chapter I] describes the different signals involved in the maritime communications 
channel, the atmospheric effects that affect them, and also demonstrates the application 
of the theory of Rician fading statistics. 

Chapter III describes the simulation program, explains the purpose of each of its 
components and illustrates some results obtained by running several simulations for 


different values of frequency and satellite’s velocity. 


Chapter IV describes the throughput of unslotted Aloha with power capture prob- 
abilitv in a Rician fading environment. 
Finally, Appendix A contains a listing of the computer program of the simulator, 


and Appendix B shows the derivation of the formulas used in chapter IV. 


Il. PHYSICS OF MARITIME MOBILE SATELLITE PROPAGATION 


A. INTRODUCTION 


Any signal in maritime mobile communications 1s typically dominated by a multi- 
path fading process caused by reflections coming from the surface of the sea. These 
multipath signals are expected to be stronger at low satellite elevation angles, where the 
main lobe of the antenna of the ship earth station (SES) picks them up. No specular 
component is expected, except for the verv calm, mirror-like sea; therefore, a model with 


Rician statistics is suggested in this chapter. 


B. MULTIPATH FADING AMPLITUDE 

It is well known that multipath reflections from the sea surface consist of specular 
and diffuse components. The relative strength of these components 1s due to wave height 
and to elevation angle. with the former the dominant factor [Ref.2]. For a smooth sea 
State the specular component dominates. As wave height increases, the diffuse compo- 
nent also increases, and over a wave height of 0.5 m, the diffuse component becomes 
constant. Sea state probability studies have been done in the past and it was concluded 
that wave heights are usuallv in a 0.5 to 5 m. range [Ref. 3]. Based on these observations. 
it is assumed that sea conditions will be rough most of the time and thus amplitude 


fading can be expressed in statistical terms. 


C. STATISTICAL MODEL FOR MULTIPATH FADING AMPLITUDE 


For a Satellite maritime link, the dominant multipath component is the diffuse 
component. [he signal can therefore be represented by two components: direct and dif- 
fuse (see Figure 1). 


O . 
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Figure 1. Propagation model for miaritinie mobile satellite contmunications 


The direct component is the line of sight (LOS) signal from the satellite. The diffuse 
component is the resultant of many random scattered waves arriving in a uniform an- 
gular distribution (Ref. 4]. If each of the waves is represented by a random phasor, the 


total difluse component at the output of the antenna 1s 


Ray= 1 explj0] = > A,expLibi] (2.1) 
i=! 
Where 
e ris the amplitude of the difluse component 


e is the phase of the diffuse component with respect to the direct component 


e <A. is the amplitude of the / scattered wave with respect to the direct component. 


If the scattered waves are sufficiently random and their phase is uniformly distrib- 
uted over 2x, Beckmann shows that the diffuse component will be Ravleigh distributed 


with a density function given by [Ref. 5] 


2 
por) = exp] - + | r20 (2.2) 


p(r) = 0 Peay 


where 
e ris the received diffuse signal level (voltage) 


© a= E{r*} 


The Ravleigh distribution requires the scattered waves to have a uniform phase dis- 
tribution over an interval of 2z. If the terrain does not provide for this requirement, the 
Ravleigh distribution will not be valid [Ref. 6]. 

The vector sum of a direct signal, with a Rayleigh fading component, results in a 
composite signal with Rician envelope statistics [Ref. 7], and the probability density 


function 1s: 





ay, r? 4 re 2Aor . 
ee (a) (2.3) 
p(r) = 90 r<O0 
ere 
® ryrepresents the envelope of the Rayleigh fading signal 
© eo = Efr} 
e A, 1s the amplitude of the direct signal 
e /,1s the modified Bessel function of zero order 
By normalizing the direct component to unity, this reduces to: 
2 
p(R) = 2B exp{ -- SET (28 R>0 (2.4) 


p(R) = 0 R<0 


where 


e k is equal to the power in the multipath component over the power in the direct 
component, that 1s: 


i ae 
2 
Ao 
¢ Ris the ratio of the amplitude of the fading signal to the amplitude of the direct 
component 
jee 
Ag 


The parameter & 1s called the Rice factor. It 1s considered a key parameter in the 
Rician distribution, and its decibel equivalent is K = 10 log, A. Figure 2 illustrates se- 


veral Rician distributions for different values of A). 


The phase of the received signal is not uniformly distributed [Ref. 5] and its density 


function is given by: 


ae - 
r(6) = = exp( -—2 Je + Gun exp(G)ltef{G)} 0<0<22 (2.5) 


where 


e Gis defined as: 
Ap cos @ 
OF 
af OF 
e erf(G) is the error function defined as: 


1 fe y 
erf(G) = i | exp| 9 Jo 
0 


Ny 27 
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Figure 2. Rician density function for ¢ = 1,f = 869 MHz, V = 0 Kts. 


1. Cumulative Distribution Function (CDF) 

The cumulative distribution function (CDF) refers either to the function G(A) 
for which the signal exceeds a level A, or to the function F(A) for which the signal 1s 
below a level A. The former is the one used in this thesis. 

The theoretical expression for the Rician CDF is obtained by integrating 


equation (2.3): 


G(A) = P(r>A) = | wre (2.6) 
A 


Figure 3 illustrates a Rician CDF. 
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Figure 3. Cumulative Distribution Function (CDF) 


Usually the CDF is plotted in Ravleigh paper which has the especial character- 
istic that a Rayleigh curve drawn in it is a straight line. This is useful for evaluating any 
outcome by visually comparing it with the Rayleigh line. Figure 4 illustrates the same 


Rician CDF seen in Figure 3 but on probability paper. 
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Figure 4. Cumulative Distribution Function (CDF) on probability paper. 


D. ATMOSPHERIC EFFECTS ON THE FADING SIGNAL COMPONENTS 


The direct and the diffuse components are both affected by Ionospheric and 
Troposphernic effects. These effects include Faraday rotation, group delay, absorption, 
dispersion, refraction and scintillation, all of which result from interaction with the 
earth’s magnetic field and the ambient electron content as the signal passes through the 
Ionosphere. Table 1 [Ref. 8] summarizes some typical maximun values obtained for an 
elevation angle of 30 °, one Way propagation, at two different frequencies. Faraday ro- 


tation can be minimized using circular polarization. Scintillation can be important at low 


latitudes (9° N to 21 ° S) and auroral latitudes (above 59 ° N), but they are insignificant 


for elevation angles greather than 10 ° and frequencies below 10 Ghz [Ref. 9]. 


Table 1. ESTIMATED MAXIMUM IONOSPHERIC EFFECTS [REF. 8]: For 
an elevation angle of 30 ° , one way propagation, and a zenith electron 
column of 10'8 electrons/? 


Frequency 
Effect 
— 850 MHz 1600 MHz 


Variation in direction of arrival 
Refraction 
Laser ane [ong Sora 










Tropospheric effects result from moisture in the atmosphere and tend to become 
more severe as frequency increases. Table 2 [Ref. 10] summarizes some estimated values 
for an elevation angle of 30 °, one wavy propagation, at two different frequencies. 

The antenna pattern will influence the magnitude of the diffuse component. The 
antenna gain rolls off below the horizon; therefore, the contribution of the diffuse com- 
ponent is due to scattered waves arriving from angles above the horizon. 

The current treatment assumes that the diffuse component is § to 14 dB (rms) below 
the direct component [Ref. 9]. 


Table 2. ESTIMATED TROPOSPHERIC ATTENUATION [REF. 10]: For an 


elevation angle of 30 °, and one way propagation 


Magnitude (dB) 
Effect - 
850 MHz 1600 MHz 


[Clear air absorption CdSSSC~iSC 
Senay) SOSC~=<“~StsOSSCSC*“‘C‘ONSC*“‘CON 
igii(moist) | 06 OT 
| Cloud attenuation | 
T Fog-attenuation ———O~dYSS 

0.05 g np (average). 0 to 73 Kmt ht. eee 

0.05 g'n (heavy), 0 to 150 mt ht. Se eee ae 
| Rain attenuation | 













Hil. SOFTWARE FADE SIMULATOR FOR A MARITIME MOBILE 
COMMUNICATION CHANNEL 


A. INTRODUCTION 


Arredondo et al. [Ref. 11] describe a hardware simulator to simulate multipath fad- 
ing for mobile radio, and Karim [Ref. 12] describes a different solution to implement this 
simulator using software. This thesis follows both reference models, adding a constant 
signal (direct component) to obtain Rician statistics. It emphasises the doppler effect 
caused by the high velocity of low altitude satellites and the vehicle’s velocity. Figure 5 


is a block diagram of the suggested model. 


B. SIMULATOR DESCRIPTION 


1. Rayleigh Multipath Fading Generator. 

The most important part of the simulator is the Rayleigh multipath fading com- 
ponent. Its implementation follows the model of Arredondo et al. [Ref. 11] and Karim 
[Ref. 12]. Two independent Gaussian sources with zero mean and equal variance are 
generated. They are passed through a shaping filter and finally added in quadrature to 
obtain a signal with an envelope that is Rayleigh distributed in amplitude and uniformly 
distributed in phase. 

The generator is illustrated in Figure 6. A key aspect in the design of this gen- 


erator is the shaping filter. 
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Figure 5. Rician fading simulator block diagram 
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Figure 6. Rayleigh fading generator block diagram 


Shaping Filter: The spectral density S(w) of the envelope of the Rayleigh fading 
signal received by an omnidirectional antenna is given by [Ref. 12]: 


4 W—- oO - 
S(w) = 3- Ul - (a se lo — Wy, < Wg (ap 


where 


* ww, 1S the carrier frequency. 


® q, is the doppler shift due to the velocity of the vehicle and the velocity of the 
satellite. 


e o? is the variance of the two independent Gaussian sequences. 


Refering to equation (3.1), taking w relative to w, and noting that we are only 
interested in relative values of the signal, the transfer function H(w) that has to be ob- 


tained to shape the input is given by: 
Way ar (3.2) 


This function is unbounded for w = @,, but following the model of Arredondo 
et al. {Ref. 11], the actual transfer function has a maximum at 0.95w, and falls to 0 with 
a slope of -18 dB / octave. It is done this wav so that the output envelope of the simu- 
lator closely fits a Rayleigh distribution. 

This transfer function will be implemented as the transfer function H(z) of a 
non-recursive digital filter with sampling frequency w, = 2q@,, and the procedure is dis- 


cussed in Karim model [Ref.12]. 


20 
H(z) =) Wanye (3.3) 
= 


where 

e IY, is the Hamming window. 

¢ a, is the Fourier coefficient of the Fourier series expansion of equation (3.2). 
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Sez is expyja@/) where f= 


Figure 7 illustrates both the theoretical and actual transfer functions. 
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Figure 7. a) Spectral density of complex envelope, b) Shaping filter freq. 


response: [Ref. 11] 


Doppler Shift [Ref 13] 

Doppler shift 1s caused by a length change with time in the propagation path. 
Let the signal from the satellite be A, cos(2nf,t). This signal arrives to the receiver after 
a time delay g(/), where g is a function describing the time variation in path length. 
Therefore the received signal is given by A,cos[2nf(tz#e(t))]. Differentiating the phase 


We can obtain the frequency of the received wave, that is: 





dg) (3.4) 


5— © [2nf(rte(s))] = ft dt 


Thus the shift frequency ts: 


dg(t) . 
Af = +fo Ti ) 








D(t) 
Since g(t) = C where 
e P(t) is the slant range, and 


e Cis the speed of light, 


we have <p i a el 

where v is the relative velocity of the receiver with respect to the source. Due to the 
higher velocity of the satellite compared with the velocity of the receiver on earth, the 
latter becomes negligible. For the purpose of this thesis, only the radial component of 
the velocity of the satellite y, will be considered. 


Therefore equation (3.5) becomes: 


ee ee (3.52) 
For satellites 
d(t) = (H+ R,) + R? — 2R,(H + R,) cos 6] (3.6) 


where 
e H is the altitude of the satellite orbit 
e 6, 1s the radius of earth = 6378.155 Km 


e 61s a linear function of time for circular orbits. 


Due to the fact that earth terminals are not located at the center of the earth, 
there 1s a maximum Doppler shift at maximum slant range (satellite “rise” and “set”), as 


shown in Figure 8. From the same figure the following are obtained: 





R,=R, +H (FF) 
R 
S21 e 4 
y = sin ( R. (3.8) 
v.=v siny (3.9) 


Where 
e 1s the total radius of the satellite orbit. 


e yis the angle described in Figure 8 


e vis the radial velocity of the satellite, from: which the maximum value of Doppler 


shift is obtained. 





Figure 8. Doppler shift froin a low altitude satellite 


Therefore, the maximum value of Doppler shift is: 


Vromax 


famax = £Sol a ) = + 7 





where 


e f,is the doppler frequency 


(3.10) 


e / is the wavelength of the carrier frequency 


e (fis the transmitted carrier frequency. 


For a circular orbit satellite, the velocity is given by: 








Oo o— ie 
Y gl Pe (Sold) 


where 


Kr 
SEG. 





COB EMeaa, + Mia) ~ BM) = 3.986013x10° 


From equations (3.9), (3.10) and (3.11) we obtain: 





Le p 
er eee 
d,max 0 CR, _ 0 CR, ee 





13.42f, 


- (3.12) 
R2 


2. Direct Component part. 
Let the direct component of the signal transmitted by the satellite be 
A, cos(w ft). The signal at the receiver will be affected by the velocity of the satellite, the 
velocity of the vehicle, or both. causing a Doppler shift. Davarian [Ref. 14] suggests a 
hardware solution to take care of the Doppler shift in the direct component: this is 
shown in Figure 10. As a result, the signal at the receiver will be A,cos{(w, + w,)t} at 
satellite rise and A,cos[(w, — w,)t] at satellite set. One of these signals is added to thie 


Rayleigh fading signal to obtain the Rician fading Signal. 
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Figure 9. Block diagram for frequency shifting of the direct path [Ref. 14 |: q, 1s 


the angular frequency of the carrier and w, denotes the Doppler shift. 


C. COMPUTER PROGRAM 
A computer prograin was developed to simulate Rician fading statistics for a mari- 


time mobile satellite conmunications channel. The computer program is written in 


Fortran and is included as Appendix A. 
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D. ANALYSIS OF RESULTS 


The simulation of several sequences of multipath Rician fading was done assuming 
a communication mobile maritime link with a satellite in circular orbit at an altitude of 
500 nautical miles. As was stated earlier in this chapter, the speed of the maritime vehicle 
is considered negligible in comparison to the speed of the satellite. As a result, the vehi- 
cle’s speed is not taken into account. Figures 10 through 12 show the generated Rician 
sequences for values of K (Rician factor) equal to 0, 5 and 10 dB, a frequency of 869 
MHz and a velocity of 30 knots. Figures 13 through 15 illustrate the histogram of the 
preceeding sequences. Figures 16 through 31 show the cumulative distribution function 
(CDF), using probability paper, of several multipath sequences for the different combi- 
nations of the following parameters K = -26, -10, 0, 10 dB, frequency = $869, 1501 
MHz, and velocity = 0, 14357 Kts. 

In order to analize the generated Rician sequences we look at two different CDF 
graphs. Comparing Figures 16 and 28 using the same frequency and velocity for values 
of K = -26 dB and K = 10 dB, it can be noticed that for the former value of K the re- 
ceived signal is greater than the abscissa for a longer period of time. This means that for 
a larger Rice factor (K) we have more fading and therefore less strength in the received 
signal at the antenna. 

Additionally, the difference in the CDF can be noticed between the curves with 0 
velocity (0 doppler) and the curves with a velocity of 14357 Kts (max. doppler); as was 


expected, the difference is larger at higher frequencies. 
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Figure 10.  Rician fading sequence for K = 0dB,f = 869 MHz, V = 30 Kts. 
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K = 5 dB, f = 869 MHz, V, = 30 Kts 
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Figure 11.  Rician fading sequence for K = 5 dB, f = 869 MHz, V = 30 Kts. 
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Figure 12.  Rician fading sequence for K = 10 dB, f = 869 MHz, V = 30 Kts. 
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K = O dB, f = 869 MHz, V, = 30 Kts 
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Figure 13. Histogram of a Rician fading sequence for K = 0dB,f = 869 MHz, V 
= 30 Kts. 
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Figure 14. Histogram of a Rician fading sequence for K = 5 dB, f = 869 MHz, V 
= 30 Kts. 


26 


K = 10 dB, f = 869 MHz, V. = 30 Kis 
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Figure 15. Histogram of a Rician fading sequence for K = 10 dB, f = 869 MHz, 
V = 30 Kts. 
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Figure 16. CDF of a Rician fading sequence for K = -26 dB, f = 869 MHz, V = 
0 Kts. 
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K = —26 dB, f = 869 MHz, V. = 14357 Kis 
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Figure 17. CDF of a Rician fading sequence for K = -26 dB, f = 869 MHz, V = 
14357 Kts. 
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K = —26 dB, f = 1501 MHz, V, = 0 Kits 
CDF 
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Figure 18. CDF of a Rician fading sequence for K = -26 dB, f = 1501 MHz, V 
= 0 Kts. 
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K = —26 dB, f = 1501 MHz, V, = 14357 Kis 
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Figure 19. CDF of a Rician fading sequence for K = -26 dB, f = 1501 MHz, V 
= 14357 Kts. 


3] 


K = —10 dB, f = 869 MHz, V 
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Figure 20. CDF of a Rician fading sequence for K = -10 dB, f = 869 MHz, V = 
0 kts. 
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K = -—10 dB, f = 869 MHz, V, = 14357 Kis 
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Figure 21. CDF of a Rician fading sequence for K = -10 dB, f = 869 MHz, V = 
14357 Kts. 
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Figure 22. CDF of a Rician fading sequence for K = -10 dB, f = 1501 MHz, V 
= 0 Kts. 
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K = -—10 dB, f = 1501 MHz, V 14357 Kts 
CDF 
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Figure 23. CDF of a Rician fading sequence for K = -10 dB, f = 1501 MHz, V 
= 14357 Kts. 
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Figure 24. CDF of a Rician fading sequence for K = 0 dB, f = 869 MHz, V = 0 
Kts. 
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K = O cB, f = 869 MHz, V, = 14357 Kis 


r 

ej0lg 
GO.G [renvcccccscececctescsceccccccsseccenscenscosscsscschsnsasasscsssusucecisasancucucoseseuetessasassuscsnsacaoscsuoncauaceugneisesecevecencncssacessscesacaee aneegencessseueeensersase 
< §g9 
WY) 
“ 
2 95 
2 Ef | nnesessnescnsoc cnsesesssosessosscoassscasosesscacosiensesecssosessees Sesssesasasssesenstossovessssensesnespacssesannessanscsiessssessosengf sacasssesesasasesesazenassceesserscsedts 
< 
E 9% 
Is, ae 
& 
nH 25 PeReUeUiaduscssossaceuecsteccecdse Masceneccesesscens 
= 10 
= S 
tJ 
ta 
~ 1 
tJ 
= Imi eae ee eee seer ce cete ed eae rene erste Gieevaccsesecracsectccesesaceceeoceort sasacnsesesaseseedonsescecsesnssecee ocsavsescesseseesieacesersseseses seats 
— 
& 0.01 en a eee occa se ccc acest os cezetescessccocscasdcccccscevacectecceessccecussoscerostccccesscccoscesecsseccescsccesesete 
30-001 a ee ee ee ee re ee Pe reate soca ccaect Jace vevesscc cues dncese cuacesecessteceocecceeesssess Ri ecevesvccsscocesetesesocvensecsessoess 


PER 





=—15 0 =o 0 =) 
RECEIVED SIGNAL RELATIVE TO LINE OF SIGHT (dB) 


Figure 25. CDF of a Rician fading sequence for K = 0 dB, f = 869 MHz, V = 
14357 Kts. 
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Figure 26. CDF of a Rician fading sequence for K = 0dB,f = 1501 MHz, V = 
0 Kts. 
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K = O dB, f = 1501 MHz, V. = 14357 Kts 
CDF 


r 


99.9 ieee: camming cecrnseececeeecensseet ecccccccccegccccecons : SQCOR GOSS ORB OODNNOOS : SQDOSONOODODOOOOOONNE eo cee cecceccccocccncs ; eovccceccccccne manne oon ccccecocccccccee: 
= QQ Jacveseereecderrececeerccenrenecertencesssrcseessneccecrssseseenassessecencccesesteeneasesseesecenaessuensaecesesaeseneaegeceeneseuscesesecaaes seaeeuenensesunensececnaacsesenesecs sees 
" 
nr QE Jorrscreredeccecscncescrscseceebecsovsessecsonseseredseeaseasersgcsscaessesecsensensen ses sseasesvenasensessosssssnesegeesstecensesseessassebaacesenseesesssesseseeasaseesenseasses 
< 90 
= 75 
aa 
- 
© 50 Melsieuislelasienslsielsieeie clele sleieineleleie'e com cleeisielsiecelc elects ceciete ceelececelecsleslecisbeesticcetcececetcecesseccesctecesciecercecses= gh cclovetesoreseveccscovessocscesocccoscsecostsssocscosovesoscecscsee 
” 

Q 25 sococoosactic 

> 

rm 10 

Ae 5 

jo 

a ee eee cance (cui cu daicesesMceauese Uceaseced: cocussincsctcdscaseseeaccsccceccctcdccsnnccctscccecescctesedcecees scgcescnsssscesstsccccscccncessostaccsesoenosceceseess 
on 1 

_ 

i 0.1 SERaEatDesicnctinescsssaemesnsecstsccat sarees rsceceeetinseeecececeececeeteeetaciscciecisccerecesiceesnsleceneccecciscccercerectionececicescrcccecececediecreccccessccecsnsepeccsenccsscscuetceceion 
<= 

o 0.01 seedeeeees maineas pe nae SQOCOCIOOCSOOOCHOONC Ls aaa SOO OCOCOCOOOOOCOCOOOC SOODDOSSOOSSN0SN0000" eecueneneccceccecseee wo ecccccccccccecoccs 
bad e ° r 4 e e e e 


é 
: 


a em ee ee ee eo a Oe ee ey ee Cee aM cee we oe We Ee gn 
SOOO SHOES OS SHES SE SHEHSSSHESSSSEH SSS HS SEHHSSHOHSHSHS HSS SOSH SS SHESSSSSHSSSHSSESHHSHSE SE SSHSSHESSSSSS SOSH SHUSESHSF SESS SHS SHSHSHH SHB SSSESCSHSSHS ES SHHS SEH SSEESE SEEK SESS SH CHK SH OHSS SHES eeFesesseeseseos 
e e e 





=o) ao °) 5 
RECEIVED SIGNAL RELATIVE TO LINE OF SIGHT (dB) 


Figure 27. CDF of a Rician fading sequence for K = 0dB,f = 1501 MHz, V = 
14357 Kts. 
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Figure 28. CDF of a Rician fading sequence for K = 10 dB, f = 869 MHz, V = 
0 Kts. 
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Figure 29. CDF of a Rician fading sequence for K = 10 dB, f = 869 MHz, V = 
14357 Kts. 


4] 


K = 10 dB, f = 1501 MHz, V 
CDF 


= O Kfts 


r 


99.9 


PERCENT TIME RECEIVED SIGNAL > ABSCISSA 
oO 
a wn 





= © 0 5 10 15 
RECEIVED SIGNAL RELATIVE TO LINE OF SIGHT (dB) 


Figure 30. CDF of a Rician fading sequence for K = 10 dB, f = 1501 MHz, V = 
0 Kts. 


x 
l 
Ss 
Q 
x 
—th 


1501 MHz, V. = 14357 Kits 
CDF 
$9.9 


99 


<5) 
90 


75 


0.1 


0.01 
0.001 


PERCENT TIME RECEIVED SIGNAL > ABSCISSA 





= 0 5 10 15 
RECEIVED SIGNAL RELATIVE TO LINE OF SIGHT (dB) 


Figure 31. CDF of a Rician fading sequence for K = 10dB,f = 1501 MHz, V = 
14357 Kts. 
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IV. CAPACITY OF UNSLOTTED ALOHA MARITIME MOBILE 
SATELLITE SYSTEMS 


A. INTRODUCTION 
Aloha 1s a multiple access method commonly used in data communications. This 
chapter analyzes the throughput of unslotted Aloha with power capture in a Rician 


fading environment. 


B. UNSLOTTED ALOHA THROUGHPUT 

In the Aloha system, a Satellite channel with a capacity of R bits per second is 
shared by a large population of users. Each user’s station randomly transmits packets 
at the channel bit rate R whenever its buffer contains one packet. When packets from 
different stations overlap (packet collision), the transmission error can be detected at the 
transmitting stations on the downlink. The stations then retransmit the packets until 
thev are free from overlap. To avoid repeated collisions, the interval of packet trans- 
mussion 1s randonuzed for each station [Ref. 15]. 

In unslotted Aloha there are two tvpes of capture effect: created and natural. Cre- 
ated capture effects occur when random power levels are used. In this case, the packet 
with the higher power level may capture the receiver when its signal-to-interference ratio 
exceeds a preset threshold. Natural capture effects occur in mobile radio communi- 
cations when varving distances between mobile nodes and base stations, as well as 
channel fading, cause varying signal powers at the base station’s receiver [Ref. 16 ]. 

As stated in the paper by Borchardt et al. (Ref. 17] about unslotted Aloha, let us 
assume an infinite number of uncoordinated users communicating to a base station us- 
ing fixed length packets. The channel traffic is approximated by a Poisson proccess with 
parameter G packets per packet length. Thus the probability that a tagged packet 1s 
overlapped by n other packets during its transmission is the arrival probability of n 
packets during the interval (4, -— 7,%+ 7) , where J is the transmission instant of the 


tagged packet and TJ is the packet transmission time. This arrival probability is 





Pr{n} = ss 


n! 


exp( —2G) (4.1) 


+n 


According to Borchardt et al. [Ref. 17], consider a receiver that captures a tagged 
packet in the presence of »# interfering packets, if the tagged packet power sufficiently 
exceeds the joint interference power. These interferers form a stochastic proccess which 
consists of early interferers, whose transmission began in the interval (t — 1, %), and late 
Mareterers witich enter the Systenmin the interval (7,7, + 1). With » interferers there are 
2” realizations of early and late interferers. C(m) 1s the number of realizations in which 
the maximum number of interferers does not exceed N (the maximum number of 
interferers that can be tolerated at any instant of time during the tagged packet trans- 


mission) and is given by [Ref. 18] 
N 
Cin)= Cn), [x] = integer>x (4.2) 


where 


neo = 7 1) 
(n) = -_s- salad 


0 , otherwise 


where Cz) 1s the number of realizations in which the maximum number of interferers 
equals /. 

The capture effect in a random access channel occurs naturally when fading reduces 
the joint power of interfering packets, as reported in {Ref. 19] for a slotted Aloha channel 
ina uniform Rayleigh fading environment. This thesis extends these results to unslotted 
Aloha in a Rician fading environment. 

It is assumed that a tagged packet will capture the receiver if its signal-to- 
interference ratio exceeds a threshold y,. 


The probability density function of a packet’s power in a Rician fading channel is 


Pas 2 
2 Aon) P Poe 
Sp(p) =k —o— ) exp( - 45 ) (4.3) 


The same expression for a normalized mean power a = | is 


folp) = exp( —AG)Io(2Agy’P ) exp( —p) (4.4) 
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Thus given ” interfering packets during the transmission of a tagged packet, the 


density function of the power of ~ interfering packets is [Ref. 17] 


n 


>» Olntfpoyy® (4.5) 


: N 
i ae 


l 
if 





Sy¥y\n) = 
Substituting equation (4.4) into equation (4.5) we get 
nr 
Avil) =se 2, Glo expl — Ap) olay) exp( —y)¥/* (4.6) 
jl] 


where j® denotes the j-fold convolution. 


Appendix B illustrates the derivation of 


| exp( — Ap)io(2Aov") exp(-y) #© for j= 1,2,3,4,5 and 6 


The density function of the power of the signal-to-interference ratio Z =< is Ob- 
tained by 
folein) = | “xfevediyerbnay 47) 
0 
Thus 


rn 
] a 2 | eee 
Folch) = LE | Wy exp(— Addio(2dov 2) exp =) D, Gol exp — Ag Io(2ov") 
- 0 j=[n/2]} (4.8) 
-exp( —y)}/*} 
The capture probability is given by 


Pr{capture|n, n>1} = 1 — P{0<Z<yoln, n=1} (4.9) 


Therefore we have 


: 
Pr{capture|n, n=1} = 1 -| “fozln)dz (4.10) 
0 
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Appendix B illustrates the derivation of the probabilities of capture for n = 1,2,3,4,5 
and 6. 
Channel throughput S is defined by 


S = GPr{capture} = G). Pr{capture|n} Pr{n} (4.11) 


n=0 


Replacing equations (4.1) and (4.9) in equation (4.11), the throughput of unslotted 


Aloha in a Rician fading channel is given by 





nN 


S= Gy 26) exp( —2G)Pr{capture|n} (4.12) 
n=Q 


and given that Pr{capture|0} = 1 1t reduces to 


GO 


Gexp(—2G)} | + by ae Pr{capture|n} (4.13) 


pr! 





i=) 


AS 


Figures 32, 33 and 34 illustrate the throughput for several capture thresholds y, us- 
ing 4) = 9.1, 4) = 0.5 and A, = | respectively, Finally, Figure 35 illustrates the through- 
put for the same capture thresholds 3. but letting parameter A, = 0 (Ravleigh fading). 
As expected the same throughput of unslotted Aloha in a Ravleigh fading environment 
eine? ] were obtained. 
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Figure 32. Unslotted Aloha throughput (A, = 0.1) 
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Figure 33. Unslotted Aloha throughput (A, = 0.5) 
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Figure 34. Unslotted Aloha throughput (A, = 1) 
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Figure 35. Unslotted Aloha throughput (A, = 0) 
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V. CONCLUSIONS 


This thesis can be divided into two main parts. The first part describes the imple- 
mentation of a software simulator for multipath fading in a maritime mobile satellite 
channel. The second part analyzes the throughput for an unslotted Aloha channel in a 
Rician fading environment. 

The fading simulator generates Rician fading sequences for any value of frequency, 
velocity and Rice Factor (K). Fades can be simulated for communication links using a 
geostationary satellite or low altitude satellites. 

Unfortunately, it was not possible to obtain data from experiments done in several 
institutions like JPL and CRC in order to compare results and make necessary adjust- 
ments to the simulator. 

For the second part of the thesis, a set of equations are derived to compute the 
throughput of unslotted Aloha in a Rician fading channel. The obtained equations are 
tested for several values of the parameters A, and y,. Finally, a comparison is made with 
results obtained for a Rayleigh fading channel [Ref. 17] by setting the parameter A, equal 
to zero, and simular results are obtained. It is noted that the throughput of a Rician 
fading Aloha channel depends on the direct component power A?. As A? increases, the 


joimt interference power also increases and consequently the throughput is reduced. 


By 


APPENDIX A. SOFTWARE RICIAN SIMULATOR 


This program generates a Rician fading sequence given variance, frequency and ve- 


locity. 

C verre §~=VARIABLE DEFINITIONS = *veevre 
C A = LOWER LIMIT 

C B = UPPER LIMIT 

C FCT = FUNCTION 

C S = PARTIAL SUM VALUE 

C H = INTERVAL LENGTH 

C AC = SIMPSON’ S RULE VALUE 

C N = NUMBER OF EVEN SUBINTERVAL 
ec NN2,NN,1I,CC = COUNTERS 

C el, = COMPONENTS OF THE RULES 
C VARIABLE DECLARATIONS 


COMPLEX RF(0: 1024) ,R(0: 1024) ,DSV(0: 1024) 
REAL SI(0: 1024) ,SQ(0: 1024) ,V,ME,G(0: 500) ,F(0: 500) ,Y,W,CONST(0: 1000 
&) ,WO,DS 
REAL*8 FR,WS,WL,A,B,FCT,S,X,H,D,AC(21),PI,Z,RGF(0: 1000) ,AX(0: 1000) 
& , IGF(0: 1000) ,MAG(0: 1024) ,TS 
INTEGER I,NN2,N,NN,NS 
© MAIN PROGRAM 
C Nl: NUMBER OF SAMPLES OF THE GAUSSIAN RANDOM GENERATOR 
N1=500 
C N2: NUMBER OF COEFFICIENTS CALCULATED FOR THE SHAPING FILTER 
N2=21 
NP = Nl + N2 


C GAUSSIAN RANDOM GENERATOR 1 WITH VARIANCE 1, MEAN 0 
ISEED=23 
WRITE(*,**) ‘ENTER VARIANCE’ 
READ(* ,**) V 
ME=0. 0 
DO 10 I=0,N1 
CALL GAUSS (ISEED,V,ME,W) 
G(1)=W 
10 CONTINUE 
é GAUSSIAN RANDOM GENERATOR 2 WITH VARIANCE 1, MEAN 0 
DO 11 L=0,N1 
CALL GAUSS (ISEED,V,ME,Y) 
BiG) —Y 
11 CONTINUE 
WRITE(8,*) (G(I),I=0,N1) 
WRITE(4,*) (FCI),I=0,N1) 


CALCULATION OF SHAPING FILTER COEFFICIENTS 


CG) C2 C2 Crm ©) 


TO CALCULATE A DEFINITE INTEGRAL FOR A FUNCTION USING SIMPSON'S 


> 


£2 


14 


LS} 


IBS, 


20 


Za 


RULE. 


NS=2i 
PI=DATAN( 1D00)*4 
C=2. 997925E+08 
WRITE(*,*) ‘ENTER OPERATING FREQ’ 
READ(*,*)FR 
WOM= 2 * Phan 
WRITE(*,*) ‘ENTER SPEED IN KNOTS' 
READ(* , ** )WM 
IF (WM .EQ. 0) THEN 
DO 1° K=0en1 
DS = COS(WO*K) 
RF(K) = CMPLX(G(K),FCK)) 
DSV(K) = CMPLX(DS,0. 0) 
R(K) = DSV(K) + RF(K) 
MAG(K) = CABS(R(K)) 
WRITE(4,*) MAG(K) 
CONTINUE 
ELSE 
WL=C/FR 
WRITE(*,*) ‘WAVELENGTH = ‘ ,WL 
WM=2**PI*(WM*1852/3600)/WL 
WS=WM/2/PI 
WRITE(*,*) 'DOPPLER FREQ = ',WS 
B=0. 977*WM 
WS=27*WM 
TS=2"*PI/WS 
WRITE(*,*) ‘ENTER (EVEN) N FOR SIMPSON INTEGRATION’ 
READ(** ,*)N 


A=0 
H=(B-A)/N 
D=H/3D00 
WRITE(* ©) ' 2ee22*COMPUTING: PLEASE WAIT=erdesececex! 
DO 13 K=0, (NS-1)/2 
S=FCT(WL,K,WS) 
DO 14 I=1,N-1, 2 
WL=A+I*H 
S=S+4. O*FCT(WL,K,WS) 
CONTINUE 
DO 15 I=2, N-2, 2 
WL=A+I*H 
S=S+2. O*FCT(WL,K,WS) 
CONTINUE 
S=S+FCT(B,K,WS) 
AC(K+1+(NS-1) /2)=S*D 
CONTINUE 
DO 20 I=1,(NS-1)/2 
AC(I)=AC(NS+1-I) 
DO 21 J=0,(NS-1) 
AX(J)=AC(J+1) 
CONTINUE 


DO 25 N=0,NS-1 
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AX(N)=AX(N)*(0. 54-0. 46*DCOS( 2*PI*N/(NS-1))) 
CONTINUE 


CONVOLUTION OF GAUSSIAN RANDOM VECTOR 1 AND SHAPING FILTER 1 
DO 26 N=N2,N1+N2 

AX(N)=0. 0D+00 
CONTINUE 


CALL CONVOL (N1,N2,G,AX,RGF) 


CONVOLUTION OF GAUSSIAN RANDOM VECTOR 2 AND SHAPING FILTER 2 
CALL CONVOL (N1,N2,F,AX, IGF) 


APLYING DOPPLER EFFECT TO DIRECT SIGNAL 
DO 30 K = 0O,NP 
DS = COS((WO + WS)*K*TS) 


RF(K) = CMPLX(RGF(K), IGF(K)) 
ADDING RAYLEIGH SIGNAL TO A CONSTANT TO OBTAIN RICEAN 
DSV(K) = CMPLX(DS,0. 0) 
R(K) = DSV(K) + RF(K) 
WRITE(8,%*) R(K) 
MAG(K) = CABS(R(K)) 
WRITE(4,%*) MAG(K) 
CONTINUE 
ENDIF 
STOP 
END 


REAL*8 FUNCTION FCTCWL,K,WS) 

REAL*8 WL,WS,PI,D 

INTEGER K 
PI=DATAN( 1D00)*4D00 
D=2D00*PI*K*WL/WS 
FCT=2D00/WS*( 1D00-( 2D00*WL/WS )****2 ) ***( -0. 25) *DCOS(D) 
RETURN 
END 


APPENDIX B. DERIVATION OF FORMULAS FOR UNSLOTTED 
ALOHA THROUGHPUT 


A. J-FOLD CONVOLUTION 


Given the expression 
2 
| eA“ o(2Aov/y Je | (B.1) 


in order to compute the j-fold convolution for} = 1, 2, 3, 4, 5 and 6, we evaluate the 


characteristic function and obtain 


a Z : 
D,{10) = | e “0 F.(2ApJSy jee dy (Bea 


—ooO 


Replacing the modified Bessel function of order 0 by its series approxirzation to the 


fourth power we Herve 


———— Te | we | 
D,{u) = ti re + e Arye” + ere 4p aan 
2 Es (B.3) 
+e” eee le 
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This 1s equivalent to 


oO 4 oO 
2 yf A a 
Ou) =e “0 | e VY ay + 43| ye dy == al yen Baap 


— ay (B.4) 
Ap ~ 3 -—(—juyy Ap 2 (i 
+ va ye “ay + 376 ye dy 


After evaluating the integrals we obtain 
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| , A A, ee 
ee he te 
(1 — ju) (1 — ju) (1 — ju) (1 — ju) 


: (B.5) 
Ao 24 
eG jx) 


Now, the j-fold convolution in time domain is equivalent to the jth power in the 


frequency domain, that 1s 
(0) P= (O (wy 


Evaluating the different jth powers and taking the inverse transform to each expression 
we obtain the desired j-fold convolutions. 
For} = 1 


4 6 


Ao Ao 
he) ~ = ae + 2Aave” + Q, 5 ee Ee + Q. 1666 Es A (B.6) 

















For] = 
(HO)? = ctl ye emer +2 =e rer +. oa ey 
8 10 12 (B.7) 
583 Ao eer Ao oy Ss 7 -y 
+ 0.5833 Sy "e + 0.1666 a aly e 
For) = 3 
4 6 
Su 3i47 ie Ao a A 
(fx)? = e “| eC. 3 e7+4.5—-yre As aye —y 
ve a ws - Ae : 
P20 yee Wel a] —y'e + 0.7083 =— yre” (B.8) 


ee ME 
+ 0.2083 aml € Wi hics 





ye + 0.0046 —— er] 


For) = 4 
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sal 1 3 = As Ap Ay 6 - 
FO)? =e sales 31 = ye 4+ +4 y" e een e* +10. 6666 = -y"e 4 
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ye” + 0.6018 ye” 
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AG An Ae = 
(Ay)? = | rT yee 05 a ee DS — ye + 20.8333 a yle? 
a> , As oa Ho) r a Ane ideal 
+ 25.8333 + 25.1666 = + 19.8611 ye 
as Pie AG 2 
+ 12.9166 ye + 6.9791 ye + 2.6157 ye” 
Ae Ane Ae 7 
+ 1.1655 pie + 0.353 ye” + 0.0848 Tyre” 
Ag? Ane Ae 
+ 0.0254 = oe 2” +.0,000128 = ory om 
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For] = 6 
- Ae _ van Z Ao ; 
(f Ay)? = aa > De) le 6 ye Lene 18 = y'e yas 36 aye y 
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eS ee feel Naot Lo 11 —-y 
45575 —— + 63.5 ye” + 61.4166 ye 
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5 5 (B.11) 
ye 20 ye Ae 
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ee ne 
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B. PROBABILITY OF CAPTURE 


First, evaluating equation (4.8), f(z |”) forn = 1 2, 3, 4, 5 and 6, the following ex- 


pressions are obtained 
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where the C,(77) values are defined in equati 
in Table 3. 
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Table 3. CONSTANT VALUES USED IN EQUATIONS (B.12) THROUGH 
(B.23) 
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Now, the probabilities of capture are evaluated by replacing equations (B.12) 
through (B.17) into equation (4.10), and we get 
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